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Density of States for Random Contractions 


Alain Joye* 

Abstract 

We define a linear functional, the DOS functional, on spaces of holomorphic functions on 
the unit disk which is associated with random ergodic contraction operators on a Hilbert 
space, in analogy with the density of state functional for random self-adjoint operators. 
The DOS functional is shown to enjoy natural integral representations on the unit circle 
and on the unit disk. For random contractions with suitable finite volume approximations, 
the DOS functional is proven to be the almost sure infinite volume limit of the trace per 
unit volume of functions of the hnite volume restrictions. Finally, in case the normalised 
counting measure of the spectrum of the finite volume restrictions converges in the infinite 
volume limit, the DOS functional is shown admit an integral representation on the disk in 
terms of the limiting measure, despite the discrepancy between the spectra of non normal 
operators and their finite volume restrictions. Moreover, the integral representation of the 
DOS functional on the unit circle is related to the Borel transform of the limiting measure. 


1 Introduction 

The density of states measure is an important mathematical notion in the study of the spec¬ 
tral properties of random self-adjoint operators, with a well dehned physical significance, 
see e.g. the textbooks [CFKSilCniKij . It is the measure associated to a positive functional 
acting on compactly supported continuous functions on the real axis, related to the random 
self-adjoint operator. Given a function, the functional is defined as the expectation of a 
diagonal matrix element of the function of the random operator, and the Riesz representa¬ 
tion theorem provides the associated density of state measure. For operators defined a 
physically appealing definition consists in considering the trace per unit volume of functions 
of the random operator restricted to boxes A C by suitable boundary conditions, and in 
taking the limit A —>• Z'^. Under ergodicity assumptions, the limit exists almost surely and 
the two notions coincide. This procedure works equally well for unitary operators, see e.g. 

HU. 

We revisit these constructions in the framework of bounded random operators that are 
not necessary normal, and study some of their properties, as described below. By rescaling, 
we can restrict attention to contraction operators that we assume are completely non unitary 
(cnu for short). 

As the there is no continuous functional calculus in this framework, we resort to the 
holomorphic functional calculus for cnu contractions developed on a Hardy space of the 
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disk, as recalled in Section [2j This allows us to define a density of state functional (DOS 
functional for short) on in analogy with that defined for self-adjoint operators, see 

Dehnition 13.11 We show in Proposition 13.21 that the DOS functional possesses a natural 
integral representation on the unit circle T by a function (p G L^(T) \ Pfg(D), and that 
when restricted to the disc algebra A(D) = U (^(I])), it further admits an integral 

representation on the disk by a complex harmonic function as Proposition 13.51 shows. 

Restrictions of random contractions to finite volume boxes A C are considered in 
Section [H under suitable ergodicity assumptions. The infinite volume limit of the trace 
per unit volume of functions of random contractions is shown to coincide with the DOS 
functional defined via the full operator in Propositions 14.31 and 14.41 We make use of this 
alternative construction of the DOS functional on A(D) to show that a priori estimates of 
the spectral radius of the finite volume restrictions imply more structure on the integral 
representations on T and B by means of p and m^f,: these functions are shown to be 
dehned by the complex conjugate of a function 'ip that is holomorphic in a neighbourhood 
of the unit disk, see Theorem ITBI and Corollary 14.81 Then, we consider the situation 
where the normalised counting measures on the spectra of finite volume restrictions of the 
random contractions converge weakly , in the infinite volume limit. Theorem 14.91 states 
that, despite the fact that finite volume restrictions of non normal operators have spectra 
that generally differ significantly from the full operator |Dll ID21 IGoKh2l ITE] , the limiting 
measure provides us with yet an alternative integral representation of the DOS functional 
on the unit disk. Moreover, the holomorphic function ijj is directly related to the Borel 
transform of the limiting measure. 

Some examples are worked out in Section [5] to illustrate the various features of the DOS 
functional. We start with random contractions defined as multiples of random unitary 
operators. Then we consider the non self-adjoint Anderson model (NSA model), whose 
hnite volume restrictions have eigenvalue distributions that give rise to limiting measure, 
see e.g. |GoKhH IGoKh2[ ID2j . Finally, the DOS functional computed for certain non 
unitary unitary band matrices, whose spectral properties are studied in m, and whose 
finite volume restrictions display similar features as those of the NSA model. 


2 Functional Calculus for CNU Contractions 

We recall from |SFBK| the main properties of the functional calculus developed for con¬ 
tractions. 

Let T be a contraction on a separable Hilbert space T-L. Consider the unique decompo¬ 
sition 

r = Tq 0 Ti on T-L = T~Lq 0 T^i, (1) 

where Tio = {p; e Ti \ \\T^ip\\ = \\ip\\ = ||r*”V’ll,^ G N}, 0 no and Tj = T\nj, 

j = 0, 1, such that Tq is unitary, and Ti is completely non unitary (cnu). The analysis of 
random unitary operators is by now well known, so we restrict attention to cnu contractions 
and therefore assume that T = Ti in the following. We recall here some basic facts from 
harmonic analysis. In the following, B denotes the open disk, its boundary is cIB that we 
will also identify with the torus T. The set of holomorphic functions on an open set S G C 
is denoted by Hol(5). 
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The Hardy class HP{0) consists in holomorphic functions on B such that 


Hull = / ®^Po<r<i if 0<p<oo, ^2) 

^ I sup^go|u(z)| if p = oo. ^ ’ 

For all 0 < p < oo, functions in HP{0) admit radial limits lim^^;^- u{re^*) = u(e**) on 9D 
for almost every t S T and In |tt(e**)| € L^(T). For 0 < p < oo, u G Ff^(B) further satisfies 
u{re^') tt(e*') in LP{T) norm. Moreover, introducing L^(T) as the set of functions in 

/ G L^(T) whose negative Fourier coefficients all vanish, the spaces Ff^(B) and L^(T) can 
be isometrically identihed, for all 1 < p < oo. These function spaces are Banach spaces, 
and even Hilbert spaces whenever p = 2. For later reference, we also introduce for all 
0 < p < oo, = {g £ HP{0) \ g{0) = 0}. Finally, the disk algebra H(B) is defined as 

the set of continuous functions on D that are holomorphic on B, i.e. H(B) = Ff°°(B)nC'(T). 
Let ~ denote the involution on holomorphic functions on B given by f{z) = f{z), and set for 
any 0 < r < 1 Ur{e'‘^) = n(re**). The following statements, among other things, are proven 
in |SFBK] . Section HI.2, Theorem 2.1. 

Theorem 2.1 Assume T is a cnu contraction on a separable Hilbert space H. Then, for 
any u G Ff°°(B), s.t. u{z) = Yln>o^riZ^ for 2 G B, u{T) is defined by the strong limit 

u{T) = s- lim UriT), where Ur{T) = '^Cnr^'T'^, (3) 


which exists. The map F1°°(B) 3 u 'o{T) G B{H) is an algebra homomorphism whieh 
further satisfies : 


a) 

b) 

c) 

d) 


u[T) = 


I if u(z) = 1 


T ifu(z) = z 

lk(7')|| < \\u\\oo 

Un{T) -3 u{T) in norm, resp. strong, resp. weak sense if 

Un ^ u uniformly on B, resp. boundedly a.e. on T, resp. boundedly on 

u{Ty = u{T*). 


( 4 ) 


Remarks 2.2 i) By the identification o/R°°(B) and L^(T), this functional calculus can 
be viewed as a homomorphism Ly^{T) 9 / —)• /(T) G with 

f(r\ = / ^ if fit) = 1 a.e. 

’ \r if f{t) = (f^ a.e. 

II/(T)|| < ll/lloo. (5) 

a) Conditions a), c) in the strong sense and m make this funetional calculus maximal and 
unique, see \SFBKf . 

Hi) If u £ ^(B), ||u„,(T) — u(r)|| —)■ 0, as n ^ oo. 

iv) If T = To © Ti, as in m, with unitary part Tq having purely absolutely continuous 
spectrum, then Theorem. \2.1\ holds with c) in the strong sense. See \SFBKf . Theorem 2.3. 
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3 DOS Functional 


We deal here with random contractions dehned on a separable Hilbert space % with 
ergodic properties we express as follows. 

We first assume some regularity assumptions. Let the probability space (H, P), where 

H is identified with {X"^ }, W C M, h G N, and P = where is a probability 

distribution on X and T is the cr-algebra generated by the cylinders. We assume that 
H 9 w !-)• G H('H) is measurable i.e. 

V(/?, V' G "H, H 9 cj !-)• (</?|T(w)'0) is measurable. (6) 

Ergodicity is expressed in the following framework. We consider % = and consider 

for j G the shift operator Sj dehned on H by 

Sjioj)k = ojk+j, k G where oj^ G X, (7) 

so that the measure P is ergodic under the set of commuting translations {Sj}j^^d. Let 
denote the canonical basis of H and Vj be the unitary operator dehned by 

Vj^k = ^k-j, VA: G Z'^. (8) 

We further assume the existence of a periodic lattice L C Z*^ spanned by {7i}i6{i,2,...,d}! 
7 i G Z*^ and the corresponding primitive cell B = Xi^i,0 < Xj < 1} n Z'^ so that for 

any j G Z'^, there exist a unique b & B and a unique <7 G L with j = h + g. 

The random contractions we consider are ergodic in the following sense: 

Ts^i^^) = VgT^V-\ V^GL. (9) 

Definition 3.1 The DOS functional L : H°°(D) C is defined for all f G by 

m = ^^Y.^{{Tb\f{TM). (10) 

I I b&B 

where \B\ denotes the cardinal of B and is a measurable cnu random contraction. 

We hrst note that 

Proposition 3.2 The map L : C is a bounded and admits the following integral 

representation: there exists ip G L^(T)\Hq(B), such that for all f G H°°(D) 

i(/) = (11) 

Remarks 3.3 0) The set L7 q(B) should be understood as the set of boundary values of these 
functions. 

i) The ergodicity assumption plays no role here. 
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ii) Such functionals on are called weakly continuous, see e.g. m Section V. 

Hi) The representation 071) says that for any G G -ffQ(]D)) 

[ /(e**)(^ + G)(e**)|^= [ V/G F-(D), (12) 

Jt Jt 

see m, Theorem 5.2. It shows that the relevant information is carried by the negative 
Fourier coefficients only. 

iv) In case ip G LP(T) C 7y^(T), with p > 1, we can actually represent L by a unique 
ip~ G obtained from (/? by substracting the contribution from the sum over positive Fourier 
coefficients. We give conditions for this to hold in Theorem \4.6\ below. 

v) Making the dependence on T of the functional L explicit in the notation, we have for all 

f G _ 

Lrif) = LT*{f). (13) 

Proof: Linearity and the bound \L{f)\ < ||/||oo stem directly from the properties of the 
functional calculus recalled above. Then one makes use of the following equivalence, see 
Theorem 5.3: L is a weakly continuous functionals on Lf°°(D) iff L is continuous under 
bounded pointwise convergence; i.e. if fn G Lf°°(B), ||/n||oo < ^ and fn{z) —)• f{z), for 
all z G D, then L{fn) — )• L{f). For such a sequence fn, we have by point c) of Thm. 12.11 
that {ipk\fn{Tu;)Fk) {ipk\f{Tu;)Fk) for all k G and \{ipk\fn{T^)Tk)\ < M, uniformly 
in n,k,uj. Hence by Lebesgue dominated convergence, we also have E(((^fe|/„(T(^)99f,)) —> 
^{{Tb\f{TLo)Tb)), for any b £ B, which yields L{fn) —)■ L{f), for \B\ finite. H 

Remarks 3.4 i) The bound \L{f)\ < ||/||oo is saturated: L{1) = 1. 

ii) For any j G N, the function z and, denoting the Fourier coefficients of 

ip by {ip{k)]k^i, 

L{J) = f = Ti-j), with (p{0) = 1. (14) 

Then we observe that for functions in 4.(1!)) C II°°(f}), we get an alternative represen¬ 
tation of L(-) on the disk. 

Let us denote by P[y?](re**) = P[(p]{x,y) the harmonic function in D given by the 
Poisson integral of y? G L^(T), with the usual abuse of notation. Due to the fact that 
lim^^;^- P[ip](re^^) = (/j(e**) almost everywhere and in L^(T) norm, we can approximate 
L{f) by an integral over smooth functions: for any / G 4(D), and any 0 < r < 1, 

L{f) = [ f{e^^)P[(p]{re^^)^+ h{r), where 

Jt ztt 

\h{r)\ < ||/||oo||</5(-)-^[¥^](r-)llLi(T) ^Oasr-^1". (15) 

Since 4(D) consists in uniformly continuous functions, we can further approximate /(e**) 
by f{re^^) = /r-(e**) to get 

L{f) = [ f{re'^)P[ip]{re"^)^ Phir), where 

Jt ztt 

Mrf < ||/||oo||¥5(-) - ^M(r-)llLi(T) + 11/-/r||H-||¥^||Li(T) ^ 0 as r ^ 1". (16) 
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The latter approximation allows us to provide L{f) with a smooth integral representa¬ 
tion over the disk. 

Proposition 3.5 Let ip G L^(T) he the integral representation of L{f) and P[ip] its Poisson 
integral on D. Then there exists a harmonic function m^p on D such that for all f G 


L{f) = / f{x + iy)mp{x,y)dxdy, where 


mp{x,y) = ^{d^ + idy){{x-iy)P[(p]{x,y)}. 


(17) 


Proof: Consider the approximation (fTUD . By Stokes theorem applied to f{z)P[ip]{z) G 

we have 

dz 


f f{re“)PM(re‘‘)L = / 

Jt jj.gn 

-I 




r‘^2i7r 

JrO OZ 

Thanks to (jl6p . we can take the limit r —)■ 1“ which yields 

m^{z) = ^^{P[ip]{z)z} = i (^[p^]{z) + (^^P[ip\{z'^ ^ . 
Using the fact that P[ip\ is harmonic, one finally gets 


(18) 


(19) 


( 20 ) 


Remark 3.6 In keeping with the fact that (p ip + G, where G G Rq(T) does not change 
the representation, one checks that 

mp+G{z) =+ G{z)/, where / f{z)G{z)dxdy = 0. (21) 

Jo 

The integral representations of L(-) discussed so far are intrinsic. There are of course 
many alternative integral representations in the disk: since A(B) C C'(B), we can extend 
L to L : (7(10)) ^ C, by means of Hahn-Banach Theorem, with ||L|| = ||L||. Given L, 
since B is compact, the Riesz Representation Theorem asserts the existence of a unique 
complex Borel measure dp, on B such that L{f) = f{x,y)dp{x,y) for all / G (7(10)) and 
||I/|| = |/r|(B). Thus, 

Lemma 3.7 Let he a measurable random contraction. There exists a complex Borel 
measure dp on B such that |/r|(B) = 1 such that 

L{f) = [_ f{x+ iy)dp{x,y), V/G A(B). (22) 

Jo 

Remark 3.8 The measure dp uniquely determined hy the extension L of L to (7(B). The 
example discussed in Section \5.4\ illustrates the fact that there may he infinitey many such 
integral representations of L. 


6 






4 Finite Volume Approximations 

Let A C be given by (2n + l)'^ symmetric translates of B along L of the form 

d 

A= IJ S + ^nj7j= IJ B + gm, (23) 

i=i ,{2n+iy 

— n<n^ <n 

with |A| = (2n + 1)‘^\B\, and let 

"Ha = span {(fj, j G A}, ='HQ (24) 

together with the corresponding orthogonal projections onto these subspaces P\,Pj^c. 


Dropping w from the notation for now, let assume that the cnu contraction T can be 
written as 

r = (25) 

where and are defined on "Ha and and is a trace class operator on "H, and 
furthermore is a cnu contraction. Such decompositions can be obtained for example by 
setting 


+ boundary conditions (26) 

+ boundary conditions, (27) 

for suitable boundary conditions at (9A, see below. For any / G 77°° (B), the operators 
/(r), f{T^) are well dehned by functional calculus and we consider two random functionals 
on 77°° (D) given by 

LkU) = ^tr(/(r^)), Za(/) = ^tr(PA/(T)PA). (28) 

From the bound ||A||i < rank(A) ||A|| on the trace norm || • ||i, we deduce that for all 
/ G 77°° (D) 

||La(/)|| < ll/lloc, I|La(/)|| < ll/lloo, (29) 

so that all arguments of the proof of Proposition 13.21 apply. Hence L\ and L\ can be 

written in the form (ttH) with corresponding random L^(T) functions <y9A(e*') and 
More precisely we have 

Lemma 4.1 Let Xj, j = 1,..., |A|, be the eigenvalues of the cnu contraction T^, repeated ac¬ 
cording to their algebraic multiplicities. Then, V / G 77°°(B), La(/) = 57 Jt f(^^*)TAiP^)dt 
where 


TA{e^') 


1 

R 


|A| 

7=1 




(30) 
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Proof: The finite dimensional contraction being cnu, (t{T\) C D. Hence, for any 

/ £ and any p <1 large enough 


LkU) = 



(31) 


For each j, since / £ Lebesgue dominated convergence implies that the limit p —)• 1 

exists, which yields the result. H 

Remarks 4.2 i) An application of Stokes theorem shows that 


LaU) 



/(^) 

(1 - xj-zr 


dxdy, 


(32) 


an expression of the fact that ^(1 — Xjz) ^ is the Bergman reproducing kernel. 

ii) Introducing the normalised counting measure of a{T^), dm^ on the closed unit disc D 

by 

dm^{x,y) =—Y^^6 {x-^X j)6{y-QXj), Aj £ (t(T^), (33) 

where the eigenvalues are repeated according to their algebraic multiplicities, we can extend 
L\ to (7(15)) by 

LA{f)= [_f{x,y)dm^{x,y), V/£ (7(10). (34) 

Jo 

However, L\{f) does not make sense for f £ (7(10)). 


4.1 Infinite Volume Limit 

Restoring the variable iv in the notation for a moment, we show that the random functionals 
LA,u]{f) and I/a,i.j(/) converge almost surely to the deterministic DOS functional L{f) as 
A —>• when / £ A(B) C By A —)• or |A| ^ oo, we mean n —oo in definition 

()23p . This is done along the same lines as in the unitary case under ergodicity assumption, 
see m, for example. 

We start by a deterministic statement: 

Proposition 4.3 Assume T and given by i25\) are cnu contractions. 

If ||PA(r — © r^'^)||i = o(|A|) as |A| —oo, then , for all f £ A(D), 


lim La(/) - 1>a(/) = 0. 
|A|^oo 


( 35 ) 





Proof: We need to show that 


^^hm^ ^ {tr(/(T^) - tr(PA/(r)PA)} = 0. (36) 

As f{T^) = P\f{T^)P\, cyclicity of the trace yields 

tr(/(r^) - Pa/(T)Pa) = tr((/(T^) - /(r))PA) = tr(PA(/(r^) - /(T))). (37) 

On the other hand, since A(B) consists in uniformly continuous functions on D, we have 
N 

P where sup |77 a^(- 3:)| — ||7?w||cx:) —y 0 as —)■ oo. (^^) 

j=o ^eS 

Hence, together with the bound ||A||i < rank(A) ||A||, we get the uniform estimate 

^|tr(i?^(r^) - PA7?iv(r)PA)| < "-^^^^WRNiT^) - Rn{T)\\ < 2||i?w||oc. (39) 

Therefore, we can focus on f{z) = z^, j G N. With 

(ri - (T^ © T^''y)p^ = ^r^(r - © r^^)(r^ © r^^'y-^-^PA (4o) 

k=0 

= ^r^(T-T'^©r'^'')PA(r'^eo)^"^"^ (4i) 

k=0 

we get for all j < N, using T and are contractions and cyclicity of the trace, 

^|tr(PA(T^' - (T^ © < j^^\\PA{T - © r^'')||i. (42) 

This estimate, the assumption ||Pa( 7~' — © T^*^)||i = o(|A|) and end the proof. H 

We finally turn to the infinite volume limit of LA,u]{f), for / G A(D), under the ergodicity 
assumption ([9]) on the way randomness enters the contraction 

Proposition 4.4 Assume is ergodic in the sense of For all f G A(]D)), 

lim TA,aj(/) = A(/), almost surely. (43) 

|A|^oo 

Proof: By construction of A, see ([251) . for any i7 : —)• C, 

(2n+iy d 

^i7(/c) = ^ ^ H{b + gm), where gm ='^ nai G T. (44) 

fceA b&B m=l i=l 
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Thus, by the ergodicity assumption, 


(2n+l)d 

La,M) = mE E i^b+gjV;j{Ts^^^^))Vg^^,+gJ 

' fceA ' fees m=i 

{2n+l)d- 

= mE E = mE E 

I I 6gB m=l I ' b&B\ni\<n 

Thanks to Birkhoff Theorem, we have on a set rij C ri of measure one, and for all b & B, 
n^JTT^ J2(^b\fiTs:;i...s:;^^u.))v’b)=n{‘fb\f{TM^ (46) 

\ni\<n 

Since ^4(0) is separable, the statement above is true for a dense countable set of functions 
{fm}meN on Cimen^fm = C n, a set of measure one. Since \B\ < oo, we infer 

T^\'^{^k\f{T^)^k) = '^T^M{v>b\f{Tu^)v’b)) a.s. (47) 

' fceA beB ' ' 


which proves the statement. 

Remarks 4.5 i) We consider A(D) only, a separable space, since H° 


is not. 




< 


— \\J Noo? 


Lehesgue dominated convergence implies that 


lim E{La^M) = , lim E 
|A|^oo |A|^oo A ^^ 

k=l 

Hi) Under the assumptions of Propositions [73| and \4U4\ 

lim LA,ui{f) = T(/), almost surely. (49) 

|A|^oo 

iv) The result also holds if T = Tq 0 Ti with purely absolutely continuous unitary part Tq. 


4.2 Spr (T^) < 1 

We show that if an a priori uniform estimate on the spectral radius of T\ holds, we deduce 
anti-analyticity of the integral representation (p of L{f). 


Let us drop the dependence on w € P in the notation. The form La is represented by 
integration against (^A(e**), see (I30p . which can be written with z = e** as 


_ 1 1^1 1 

<Pa{z) = iPa{z), where Pa(2) = E 

i= ’ 


(50) 


1^1 

As Ta is cnu, pA is holomorphic in D and we have the absolutely converging power series 


V'a(^) = E 


oo , m* i’ 

tr Tl k w 

^ z^, Vz S 


k=0 


|A| 


(51) 
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Theorem 4.6 Let satisfy the hypotheses of Propositions and \4.4\ Assume there 
exists r < 1 such that 

Spr (Ta) < r, y andy u en. (52) 

Then, there exists € Hol{9/r) such that Proposition Id.SI holds with = fj{e^^): 

m = l^f{e“)W^)L. (53) 

Remarks 4.7 i) We can satisfy the hypothesis by properly rescaling the operator T^. 
ii) The property Spr(TA) < r < 1 VA does not imply Spr(T) < 1. Indeed, finite volume 
approximations of non normal operators typically miss important parts ofa(T), see Section 
El and ^8^ in particular. 


Proof: If (1521) holds, is holomorphic in the larger disc D/r for all A G oj ^ Ll and 

\'f’Aiz)\ < -—Vz G D/r. (54) 

1 — r\z\ 

In particular, the family {V'aIa of holomorphic functions on D/r is uniformly bounded 
on each compact subset of D/r. Hence, by Montel Theorem, see e.g. [Rj, Theorem 14.6, 
{V^aIa is a normal family. Therefore, for each hxed cj G Hq) the set of measure one on which 
Proposition 14.41 holds, there exists a subsequence {V'Ai.jfceN which converges uniformly on 
each compact subsets of D/r to a function 'if{z) which is holomorphic on D/r. In particular, 
for all / G A(D), 

lim LAfc(/) = lim f /(e**)V'Afe(e**)|^ = / /(e**)?/(e*‘)|^, (55) 

k^oo k^oo Jy 7t 

where ipiz) is analytic in a neighbourhood of D. By Remark 14.51 iiii. we get that ip G L^(T) 
which represents L{f) is given by (p{t) = and if is independent of w. H 

Consequently, 

Corollary 4.8 With fj = ^ Ro/(D/r), and for all f{z) = ^ 

OO 

L{f) = (^I/)l^(t) = (^l/)p(N) = '^Kan- (56) 

n=0 


The integral representation reads 


m^{z) = m:^iz) 



(57) 


4.3 Representations of L via finite volume approximations 

Let us consider now the random counting measure dm^ (l33]) and assume that it admits a 
weak limit, almost surely: for all a; G Hq with P(Ho) = 1; and for all / G C'(D), 

hm LA,a.(/)= lim l_f{x,y)dmt{x,y)= if{x,y)dm^{x,y). (58) 
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Then, for any a; G Oq) duii^ > 0 provides another representation of L on A(D), since, 
specialising to / G ^(O), we get from (f4^ 

L{f)= [_fix + iy)dmu;{x,y). (59) 

Jd 


Theorem 4.9 Let T^, satisfy the assumptions of Propositions [7^ and \4-4\ o.nd dm^ in (3^ 
converge weakly to drui^, almost surely. Then L admits the following representation 



f{x + iy)dm{x,y), 


V/ G ^(B), 


(60) 


where dm > 0 is given by E(dm(j). 

Further assume Spr (T^) < r < 1, for all ui ^ LI and A. 


Then, Theorem 4-6 holds with 




and Proposition \3.5\ holds with 


dm{x, y) 

1 — z{x — iy) ’ 


Vz G B/r, 


If dm{x,y) 

^ 'x Jo [I - z{x + iy)y 


(61) 


(62) 


Remark 4.10 If the weak limit of the normalised counting measure of the finite volume 
spectrum exists, see e.g. jGoKhSlJ for such cases, it provides another representation of the 
DOS functional. In that sense, the spectrum of the finite volume restrictions acquire a global 
meaning, in spite of the fact that it can be very different from the spectrum of the infinite 
volume operator. In particular, the support of the limiting measure can be disjoint from the 
spectrum of the operator, see Section^^f^for such an example. 


Proof: The first statement is a consequence of (1591) and Fubini’s Theorem. The assump¬ 
tion on Spr(T^) implies supp dm C {\z\ < r}. Then, Corollary 14.81 applied to f{z) = z^ 
yields the coefficients of the power expansion of the holomorphic function 

y{z) = E bkz’^, z G D/r. (63) 

k>Q 


We have 

bk = lix — iy)^dfh{x,y), with \hk\ <r^, k £ N. (64) 

Jn 

Thus, exchanging integration and summation, we get expression (1611) for z G D/r. The last 
statement follows from Corollary 14.81 H 

Remark 4.11 Thanks to Remarks 12.21 iv). and \4-5\ iv), all results of Sectionf^hold if T 
writes as T = Tq © Ti, see m, with a unitary part Tq that is purely absolutely continuous. 
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5 Special Cases 

We take a closer look at various particular cases allowing us to get further informations on 
the integral representation ip. 


5.1 The Normal Case 

A first special case of interest occurs when is normal, i.e., when there exists orthogonal 
projection valued measures {dEi^{x + iy)}{x+iy)&a{T^) such that in the weak sense, 

Toj= [ {x + iy)dE^{x + iy). (65) 

In such a case, we have a continuous functional calculus: for any / G 

/(TL) = / f{x + iy)dE^{x + iy). (66) 

Hence, assuming that is normal for any a; € fl. Definition 13.11 gives rise to a positive 
functional on (^(ll])), so that by Riesz-Markov Theorem 

L{f) = f{x + iy)dm{x,y), where dm is a non-negative Borel measure on B. (67) 

Jo 

In this favourable framework, we have 

Lemma 5.1 LetT^^T^^ he contractions such that \\P\{T—T^(BT^^)\\i = o(|A|), uniformly 
in (jj G D. Further assume T and are normal and the ergodicity assumption (0) holds. 
Then, as A —)• oo, 

dm^ —>■ dm a.s., in the weak-* sense. (68) 

Proof: The same arguments using Stone Weierstrass and Birkhoff theorems together with 
the separability of (7(15)) prove the result as in the previous section. H 


Applying Propositions 13.21 and 13.51 to the normal case, we get for any / G A(B) 


/ f{x + iy)m^{x,y)dxdy = / f{x + iy)dm{x,y), 
JO Jo 


(69) 


where dm{x,y) is the non negative usual density of states, and m^p{x,y) is harmonic and 
in general complex valued. This special case makes explicit the lack of uniqueness in the 
representation of analytic functionals. 


Remark 5.2 If, moreover, there exists 0 < r < 1 such that for all lo G D, Spr {T^^) < r 
and Tcj is cnu, then Theorem \4.S\ holds with dm = dm, the density of states. In particular. 


Hf) = 


f{z)dm{x,y) dzAdz 
(1 — z(x I-iy))^ 2i-K ’ 


Vz G D/r. 


(70) 
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5.2 Multiple of Unitary Operators 

Consider now a special normal case where the statement above allows us to make the link 
between the DOS functional and the density of state measure of a random unitary operator 
C4; more explicit. Let 0 < r < 1 and U^j be a random unitary operator defined on with 

w G T^, which is measurable and ergodic. Models of this type are studied in [BHJlUT] . The 
details do not matter for our purpose here. We consider the random normal cnu contraction 


U = rU^. (71) 

We assume that there exist finite volume approximations such that is a finite dimen¬ 
sional unitary matrix; see e.g. m for examples of this situation with A = {—n + l,n}. 
Consequently, we have the trivial bounds Spr(T^) = Spi (Ti^) = r, r < 1. Stressing the r 
dependence in the notation, the DOS functional on A(D) is represented by 

^(r)(j) _ f where G Hol(B/r). (72) 

Jt 27r 

The density of states measure dk for ergodic unitary operators is a normalised positive reg¬ 
ular Borel measure on T characterised as in Section [5. II by Definition 13. II with / continuous 
on the circle via Riesz-Markov Theorem, m- 

2 

[ fie^')dk{t), V/gC(9B). (73) 

k=i df 

The explicit link is provided by 


Proposition 5.3 Let = rlJuj, with 0 < r < 1. With the notations and assumptions 
above, for all \z\ < Ijr, 


In particular, 


and, 




dk{t') 

1 — zre~'^^' 


(74) 

'iPre^d-t'Y 
1 — 

dk{t'), 

(75) 

]^P[dk]{rP^) > 0, 

(76) 


where P[dk]{z) denotes the Poisson integral in B of the non negative measure dk on T. 


Remarks 5.4 i) The representation ^7^ shows that coincides with the Borel (or 

Cauchy) transform of the measure dk on T, taken at point rz G B. The function may 
admit analytic extensions outside B/r, depending on dk. 

a) While Toj = is cnu, the boundary values as r —)• 1“ of the real part of the integral 
representation of L{f) yield the absolutely continuous component of the density of state 
measure ofUi^. Indeed, if dk{t) = ^(e**)^ + dp.s{t) is the Lebesgue decomposition of dk, 

lim — 1 = ^(e**), almost everywhere on T, (77) 

r—)-l~ 


see Theorem 11.24- 


14 









Proof: Let / G ^(D). From ([72|) and (1731) . we get 


[ = [ f{re^^)dk{t). (78) 

Jt Jt 

The coefficients of the expansion for 2 G D/r are given by = 

r'^k{n),n G N, where k{n) = fj dk{t)e~^^'^. Hence, with sign(O) = 0, 




int 


riGN 


r - 1 


—i Ylnez sign(n)rl"'lfc(n)e*”*. 


(79) 


The last two series coincide with the real and imaginary parts of dk * {Pr{') + iQri'))iP'^), 
where Pr{t) and Qr{i), are the Poisson and conjugate Poisson kernels given by 


Pr{t) + iQr{t) 


1 + re** 
1 — re** 


1 — r^ + i2r sin(t) 

-5—^-rT> 0 <r<l 

1 _l_ r2 _ 2r cos(t) 


(80) 


Replacing e** in (1751) by 2 ; G D/r yields d74l) . 


5.3 Non Self-Adjoint Anderson Model 

The non self-adjoint Anderson model (NSA model for short), provides an example in which 
the finite volume and infinite volume versions of the random operator have quite different 
spectra that can be computed explicitly. See |D2] and [GoKh2] for more general non self- 
adjoint random operators with similar properties. After suitable rescaling, the NSA model 
provides us with an illustration of our results. 

The NSA model is a one parameter deformation of the one dimensional random 

Anderson model of solid state physics defined as follows in the canonical basis of f^(Z): 

+ iOjifj, Vj G Z. (81) 

We assume that g > 0 and {ujj}j^z are i.i.d. real valued random variables distributed 
according to a measure dg supported on a compact interval [—B, B], Let Eg be the ellipse 

Eg = {e**+*^ + I 0 G [0, 27r]}, (82) 

which coincides with the spectrum of P[Q{g), the NSA model in absence of random potential. 
It is proven in |D2j that, provided B > , 

a{Hi_j{g)) = Eg + [—B,B\, almost surely. (83) 

Moreover, consider H^{g) the finite volume restriction of H^{g) to /^(A), where A = {—n, n} 
with Dirichlet boundary conditions. Then, for S' > 0, the matrix H^{g) is similar to H^{0) 
with similarity transform given hy W so that its spectrum is real. The matrix 
iL/)(0) is the finite volume restriction of the Anderson model 77^(0), with 

= Eq + [—B, B], almost surely. (84) 
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To cast these considerations in our framework, we assume B > + e ^ and set 

Tuj{g) = H^{g)/{e<^ + e"5 + B), resp. T^^ig) = H^{g)/{e3 + e~<^ + B), (85) 

on resp. /^(A) with Dirichlet boundary conditions. By construction, both operators 

are contractions. Moreover 

Lemma 5.5 The matrix T^[g) is cnu for |A| large enough, and all u £ Q, whereas the 
operator T^{g)) is cnu almost surely. Also, for g > 0, 

Spr (r^( 5 r)) = 1, a.s. and lim Spr {T^{g)) = ^ o < Vw G O. (86) 

Proof: Statements (15H1) are consequences of (15^ . (IMl) . and properties of hnite volume 

approximations of self adjoint operators. Thus T^{g) is cnu for |A| large enough. Given 
([5^ . Ti^{g) is cnu if ±1 are not eigenvalues of T^{g). Suppose ip± G l‘^{Z) are normalized 
and satisfy Ti,j{g)ip± = ±ip±. Then, using the fact that T*{g) is a contraction, and Cauchy- 
Schwarz inequality, we get T*{g)ip± = dsnp±. Hence ±1 are eigenvalues with eigenvectors 
ip± for the self-adjoint Anderson type operator (Ti^{g) -|- T*{g))/2. But this is known to 
happen with zero probability only. H 


Stressing the g dependence in the notation, we denote the DOS functional on A(B) by 
^(9)(.). The foregoing and Theorem 14.91 immediately show that 

Lemma 5.6 Let dm^{g) be the normalised counting measure of a{T^{g))). Then, for any 

5 > 0 , 

dm^{g) —>■ dkg, a.s., in the weak sense, (87) 

where dkg > 0 is the density of states of the rescaled self-adjoint Anderson model ■ 


Hence, for any f G A(]D)), and any g > 0, 


£(»)(/)=/■ 

Jl-2-eB,2-eB] ye^ + e 9-\-B J 


i[-2-eB,2-eB] 

where dk > 0 is the density of states of H^{0). 


dk{x), 


( 88 ) 


Remark 5.7 The support of dkg is a subset of the almost sure spectrum ofT,_g{g). 


The estimates provided in Lemma 15.51 show that L^9)^.'j admits an integral representation 
on T in term of a holomorphic function ^(^^(z), see Theorem 14.61 Moreover, 

Lemma 5.8 Set s{g) = e9 e~9 -|- B. The function 'tp^9) {g analytic and admits 
a converging power series in {\z\ < s{g)/{2-\-B)}. Moreover, 'ijj^9) admits an analytic 
continuation on C\ {\s{g)/{2 -\- B),oo[U] — oo, —s{g)/{2 -\- R)]} given by 


^(»)(.) = f 

Ji-{ 2 +b), 2 +b] 1 - zx/s{g) 


(89) 
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Remark 5.9 With the Borel transform Fk of dk given by 


Fkiz)= [ V2GC\[-2 + R,2 + R], 

Jl-2+B,2+B] X — Z 


we have for z ^ 0, 

■0(9) (z) = -Fk {s{g)/z) s{g)/z. 


(90) 


(91) 


Proof: It is a special case of Theorem 14.91 where the coefficients of the power expansion 
of are given by 


bn — bn 


I[-2+B,2+B] X'^dk{x) 

{e9 + e-9 + BY ’ 


, \bn\ < 


(2 + R)^ 
s{gY 


Vn G N. 


(92) 


5.4 Non Unitary Band Matrices 

We further illustrate both the discrepancy between limiting measure and spectrum of the 
full operator, and the multiplicity of representations of the DOS functional by measures on 
B by considering random contractions with a band matrix representation. 


Let be defined on /^(Z) by its matrix representation in the canonical basis given as 






0 

0 



0 

0 



ei‘^2j+2^ 


0 

0 



0 

0 


\ 


V 


giaj2j+4Q, 


giaj2j+4 p ' ■ J 


(93) 


where {ujj}j^z are T-valued iid random variables. The deterministic coefficients a,f3,'y,5 
are assumed to give rise to a non unitary matrix 


Co 


(“ fj s.t. ||Co||c2<l. 


(94) 


Actually, = D^^T, where is a diagonal unitary random operator with elements 
and T is a deterministic operator whose representation has the form (|93l) . where all ojj’s 
are equal to zero. Such random ergodic operators arise in the analysis of certain random 
quantum walks and some of their spectral properties are analysed in m- In particular, it 
is shown there that is a cnu contraction on /^(Z) if and only if 


detCol < 1, |a| < 1, and |(5| < 1. 


(95) 


Moreover, HUH = 1, and Spr(Ttj) may take the value 1, depending on the parameters. 
We hrst note the following simple general result: 
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Lemma 5.10 Let defined by Ii93\) with (94-), and Ii95\) . and assume the random phases 
are uniformly distributed. Then 


L(/) = /(0), V/e^(D). (96) 

Proof: One first observes that for any n € N*, any fc € Z, {ipk\T)fiipk) = 
where T{uj{k)) is a random variable independent of uj^, and m{k) G N*. This is a conse¬ 
quence of the fact that no cancellation of the random phases can occur due to the shape 
(f93]l of the matrix representation of Hence, = 0, for all A; G Z and n G N*. 

Thus, approximating any / G H(]D)) by a polynomial, we get the result. H 


Remark 5.11 We get from I196\) that the integral representations of L{-) on T and D are 
given by tp = 1 and ^ respectively. 

The DOS functional (|96p can be extended to (7(19)) as integration against a Dirac measure 
at the origin, but not only. For example, using in polar coordinates, defining for any 

0 < p < 1 

= IT f- f{r cos{e),r sm{6))6p{r)de = [ f{x,y)dpl^P\x,y), (97) 

we get a one parameter family of extensions of L corresponding to integration against 
positive measures on D, with support of the circle of radius p, centered at the origin, 
such that L)p){f) = /(O) if / G H(B). The support of dp,)^) may or may not belong to the 
spectrum of the examples treated in m show that for det Cq = 0 and uniform i.i.d. 
phases Uj, the spectrum is given by the origin and ring centered at the origin whose 
radiuses depend on the parameters. Note finally that any radial probability measure with 
smooth density p{r) provides us with an extension L(-) on C'(]D>) of the form 

L{f) = [ f{rcos{6),rsm{9))p{r)rdrd6, (98) 

Jo 

which agrees with L(-) on H(D), thanks to the mean value property of harmonic functions. 


We now specify the values of the parameters to Cq = , 0 < 5 < 1 , see (IM)l . and 

look at finite volume restriction generated by finite rank perturbations. Note that in this 
case, is not cnu since 

~ S © gS, (99) 

where ~ denotes unitary equivalence and S is the standard shift on /^(Z). Remark 14.111 
shows we can nevertheless apply our results to this case. 

Setting A = { 1 , 2 ,... , 2 n}, we dehne = D^T^, with unitary and diagonal and 

/O 0 5 \ 


rj-^A _ 


1 0 0 
0 0 

1 0 


0 9 
0 0 

0 0 g 

1 0 0 / 


( 100 ) 
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This yields a finite volume restriction of T^, that is a cnu contraction with eigenvalues 

Aj(a;) = j = l,...,2n. (101) 

The factor ^ tends to 0 as n —)• oo almost surely, by our assumption on the 

distribution of the phases, and we get for any / G ^(ID)), 

2n 

lim = / f{V9cos{e),y/gsm{6))^, (102) 

n^oo 2n ^ 7(0,27r] 27r 

by a Riemann sum argument. In other words, the normalised counting measure on ct(T^) 
converges weakly to dm = 5 in polar coordinates. As Spr(T^) = < 1, we 

compute from Theorem 14.91 

S 1 - "*■ 1^1 ■= 

SO that we recover (/j(e**) = 1, in keeping with Bemark l5.11l 

Remark 5.12 The support of the limiting measure dm is disjoint from a{Ti^) in this case: 

supp dm n (7(T^) = s/gdn n {^(SD U dD} = 0. (104) 
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